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Abstract 

An iterative procedure is developed with the aim of constructing homogeneity 
rules for the distribution P(p, 5) of the particle density p at resolution scale 5. 
A single iteration step consists of a change in the normalization point of 
P(p, 5) followed by a rescaling. Similar transformation rule is introduced 
for density fluctuations contaminated by Poisson noise. Application of the 
iterative procedure is given for the Ginzburg-Landau description of phase- 
transition from the quark-gluon plasma and for random cascading models. 
PACS numbers: 13.85.Hd, ll.10.Gh, 05.40.+j 
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Over the past ten years significant changes have occured in our knowledge of how to 
analyse the multiparticle density fluctuations and correlations, see Ref. |I] for a recent review. 
Several new techniques have been developed dominantly with the same goal: to search for 
self-similar patterns possessing no characteristic scale over a certain range of resolving power. 
Such fluctuations are expected to occur e.g. in elementary collision processes due to scale- 
invariant branching of parton cascades or in heavy-ion collisions as remnants of the phase 
transition from quark-gluon plasma to hadronic matter. 

Self-similarity manifests itself in the power-law dependence of quantities characterizing 
the strength of fluctuations or correlations as the resolution scale (bin size in phase-space) 
is varied Q. Besides power laws, there is another possible source of scaling behavior in 
multiparticle production. The probability distribution P(p, 5) of the particle density p at 
resolution scale 5 can be a generalized homogeneous function of its arguments, P(p, 5) ~ 
5 _/3 JF(m = p/5 u ) with j3 and v being certain parameters [1]. In this case scaling refers to the 
data collapsing feature of homogeneous functions: after suitable rescaling the data points 
corresponding to different distributions P(p,5) collapse onto a single curve, F{u). 

Homogeneity rules play a central role in the theory of critical phenomena. Near the crit- 
ical point of a physical system the thermodynamic functions exhibit homogeneous form |2| 
which implies the existence of a law of corresponding states. Using a suitably chosen scal- 
ing transformation it is possible to bring different states of the same system to coincidence 
and thus to compress many experimental or theoretical results into a compact form. This 
scaling hypothesis was first justified by the spin decimation procedure [||], the prototype of 
renormalization group (RG) transformations in statistical physics f|,[5|. 

The connection between scale-invariant density fluctuations in multiparticle dynamics 
and the scaling laws of statistical physics has already been established 0,0. In this Letter 
an attempt is made to deepen the relationship between the two fields by developing an 
iterative procedure that acts on the distribution P(p, 5) and provides the analogue of RG 
transformations. In the course of the iteration one may (and will) arrive at scaling behavior 
of P(p, 5) even if the distribution is not a homogeneous function of its arguments at the 
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beginning of the procedure. 

Let us start with some elementary properties of the probability distribution P(p, 5) of 
observing particle density p at resolution scale (bin size) 5. For simplicity we shall consider 
here only averaging over the events. The effect of statistical noise will be discussed later on. 
The probability density P(p,S) fulfills the normalization condition / °° P(p,S)dp = 1. Its 
moments are defined by (p q (5)) = / °° p q P(p,5)dp. If only the first moment (p(S)) is respon- 
sible for the £- dependence of P(p, 5) the distribution can be written in the homogeneous 
form 

where the function T depends only on the scaling variable u = p/(p(5)). It satisfies the 
two normalization conditions / °° Fiujdu = / °° uT[u)du = 1. Thus we have two properly 
normalized scaling functions, the second one is F(u) = u!F(u) = pP(u). Obviously the 
normalized moments of P(p, 6) coincide with the moments of T(u), 



roa 

/ u q F(u)du = C q (2) 
Jo 



(p(5)) q 

which are independent of 5. Moreover, the moments of T{u) are given by C q = 
Jq 00 u q iF(v)du = C q+ i i.e. the difference between the two scaling functions amounts to a 
shift in the ranks of their moments. 

Although the scaling law of type Eq. (P is frequently encountered in nature, for self- 
similar multiparticle density fluctuations it can not be valid. The observations indicate 
power-law dependence of the normalized moments of P(p, 5) as the resolution scale is varied, 
C q (5) oc where the (p q are the intermittency exponents [|]]. Accordingly, we have 

(^-dependent "scaling" functions T{u, 5) and ^F(u,5). Let us observe that the latter enables 
one to perform a rescaling of type Eq. (0) yet again, the new scaling hypothesis reads 

6) = -J— T x ( J^—) . (3) 

(Rescaling JF(m, 5) in a similar manner has no effect since C\(5) = 1.) The subscript of the 
new scaling function refers to the first step of an iterative sequence of transformations acting 



on the distribution P(p,5). For the starting distributions and moments corresponding to 
the Oth step we shall not indicate this explicitly. 

Our new scaling variable is u\ = u/C\{5) and the first-iterate moments C q> i = 
Jo°° u\ J r i(ui)dui are related to the original ones through 

C q (S) C q+ i(5) 

Expressing now the scaling function J-i{u\) similarly we get J-'iiui) = C2(5) pP \u/Cz{S)\ 
which can be rewritten as 

It is seen that J-x{ui) has the same form as J-{u), namely, J-i(ui) = (p(5))i P±[p/ {p(S))i} 
where the first-iterate distribution P\(p, 5) and its moments (p g {S))i are obtained by chang- 
ing the normalization point of P(p, 5) from the Oth moment to (p{5)) and rescaling by (p{5)) 
to maintain the overall normalization. 

Demonstrating the benefits of the above transformation of P(p, 5) let us assume that 
the self-similar density fluctuations show monofractal properties. Then the intermittency 
exponents are given by <p q = (p2(q — l) and the anomalous fractal dimensions D q = ip q / (q— 1) 
are g-independent, D q = Z^- The normalized moments of P(p,5) take the form C q (8) = 
A q [C2(8)] q ~ 1 with coefficients A q independent of 5. However, Eq. @ yields A q = Cq-i^. 
This implies that monofractal density fluctuations give rise not only to power-law scaling of 
the normalized moments of P(p, 5) but also to homogeneous form of the distribution after 
changing its normalization point to (p{5)) and making an appropriate rescaling. 

In the followings we omit explicit reference to the variable 5 to simplify the notation. If 
^i(wi) is not independent of the resolution scale, the single iteration step described before 
can be repeated until the appearance of scaling (i.e. data collapsing) behavior provided 
that certain moments of P(p) exist, see below. In the kth step of the iteration we have 
the normalization conditions f£° J-k(uk)duk = Jq° Tk{uk)duk = 1 where Tkiuk) = UkJ-k(uk)- 
The distribution Tkiuk) = (p)k Pk(p/(p)k) is related to the (k — l)th step through 



Fk{u k ) = Ci jfe _i T k -\ ( ~ k 1 I (6) 



vCi,fc-i 

and its moments C q>k = Jo°° ' ul^kiu^duk are given by 
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The (k + l)th step of the iteration follows in the same way as the first step in Eq. (^), 

Fk(u k ) = T k+1 [^-\ . (8) 

Scaling behavior of P k+ i(p) appears if the monofractal type relationship C q>k oc C^ 1 holds 
with constants of proportionality independent of the resolution scale S. 

Let us now consider how the distribution P k (p) is related to the original P(p). The 
connection is provided by the straightforward extension of the normalization point changing 
transformation discussed previously, 

and thus the moments of Pk(p) are given by {p q )k = (p q+k )/(p )■ From this we get for the 
normalized moments C qjk the expression 

C q , k = || = C q+k ^- (10) 

in terms of the original normalized moments C q . Homogeneity of P k +x(p) appears in case of 
^-independence of the product Ck +q C^2^k+v The necessar y condition of performing the 
kth iteration step is the existence of the moments of P(p) up to kth order. 

Eq. (§) is analogous to the RG transformations which correspond in many important 
cases to a change in the norm of the parameters characterizing a physical system ||. It is 
a Gell-Mann and Low type relationship [0] in which a multiplicative transformation (in our 
case multiplication by p k ) is compensated by a rescaling and a suitable change in the renor- 
malized parameters. Since the operation P k (p) — » -fV-i(p) does not exist the iteration steps 
constitute a semi-group. In statistical physics the Kadanoff- Wilson type RG transformations 



perform a systematic reduction in the number of degrees of freedom (spin decimation). This 
enables one to eliminate the small-scale fluctuations from the problem which are irrelevant 
to critical point behavior such as the homogeneity of thermodynamic functions. In the itera- 
tive procedure developed here the elimination of fluctuations manifests through elimination 
of the moments of P(p) via moment-shifting and rescaling until one observes a homogeneity 
rule of type Eq. ([!]) for Pk{p)- 

The finite number of detected particles gives rise to statistical fluctuations contaminating 
the distribution P(p) of purely dynamical origin. It is convenient to represent the statistical 
component by a Poisson process whose convolution with P(p) can be described by the 
Poisson transform |IJ. It yields the discrete probability distribution P(n) of detecting n 
particles: 

P(n) = / P{p) P -e~"dp. (11) 
Jo n\ 

The Poisson transform has the useful property that the factorial moments £ 9 of P(n) are 
equivalent to the ordinary moments of P{p), 

oo 

i q = {n^) = Y j n [q] P{n) = {p") (12) 

n=q 

where = n(n — 1) . . . (n — q + 1). Thus in the presence of Poisson noise factorial observ- 
ables make it possible to study directly the dynamical fluctuations whose self-similarity is 
signalled by the power-law dependence of the normalized factorial moments F q = on 
the resolution scale 5 0. 

We are interested here in the Poisson transform of the iterated dynamical distribution 
Pfc(p). But let us first consider briefly some properties of Eq. @ itself. This transformation 
of P(p) is known as size-biasing of order k in the mathematical literature [Kj. A worth 
noting feature is that the log-exponential family of continuous probability densities are 
form-invariant under size-biasing, only their parameters are affected. In this respect the log- 
normal distribution deserves special attention. Its parameters are changed by size-biasing 
of order k according to p — > p + ka 2 with p and a 2 being the mean and variance of the 
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original distribution. Therefore a log- normal P(p) is a fixed point of our iterative procedure: 
although Eq. (||) changes its first moment, this change is scaled out by constructing !Fk( u k) 
and we get an unaffected scaling function at each iteration step. 

The Poisson transform of Pk{p) leads to size-biasing of factorial order k acting on P(n): 

poo n n 77,^] 

W=l P k (p)^e-Pdp= J - JM -P(n) (13) 

\ ' <=k 

where the symbol -<=A; denotes left-shifting the size-biased P(n) by k places. It is seen that 
we have the same kind of transformation as for P(p): the normalization point of P(n) is 
changed to (n' fc ') and a rescaling by (n^) is performed to maintain the overall normalization. 
This yields a discrete distribution with k leading zeroes, P(n < k) = 0, hence a left-shift by k 
places is required to maintain also the original support in n of the nonzero probabilities |10| . 
The normalized factorial moments F q ^ are related to F q in the same way as Cq,k to Cq 
i.e. F q ^k = Fq+kPj!' 1 Pk+i- The expression for Pk{n) implies that the normalization point 
changing transformations form a discrete semi-group in the presence of Poisson noise. On the 
contrary, Pk{p) is well defined for fractional k too although the iterative procedure that led 
us to Eq. (^[) enables only integer valued k. It is worth mentioning that F 2) k-2 — PkPk-2Pk-i 
yields for k > 2 the kth bunching parameter |TT| in the limit <5 — > 0. 

We now turn our attention to applications of our results. In recent years considerable 
interest has been devoted to the Ginzburg-Landau theory of the phase transition from quark- 
gluon plasma to hadronic matter |12]JT3[] . Instead of a strict power-law scaling F q oc 5~' Pq 
of the normalized factorial moments, a so-called F-scaling behavior F g oc Ft is observed 
with f3 q = (q — l) u . For second-order transition v = 1.304 and in the case of a first-order 
transition v = 1.45. The F-scaling rule of the above form has the property Fq ^ 1 hence 
it is difficult to identify the family of discrete probability laws that belongs to the observed 
Fg-pattern. Performing a single iteration step the factorial moments change according to 
F g:1 oc F2 9 with (3 g = q u — q and therefore 

F^ocF-g with &=|£f. (14) 
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This form of F-scaling is the log-Levy law |15j with Levy index < v < 2. We have thus 
identified the family of probability laws for P\{n) in the Ginzburg-Landau formalism with a 
strict bound on v. The theoretically relevant values are within the allowed range. Our result 
can be checked in a laser experiment at the threshold of lasing where the photocount fluctu- 
ations are in accordance with the Ginzburg-Landau theory of second-order phase transitions 
and the F-scaling rule with v = 1.304 has been experimentally confirmed fl4j . 



The log-Levy law F q oc with (3 q given in Eq. fli"4] ) is introduced to describe multiplica- 



tive random cascading models ||15|| . Let us keep track of the evolution of factorial moments 
in the course of successive iteration steps. It can be characterized by the ratio 

D q ,k _ (q + ky + (q-l)k»-q(k + iy 

D q q"-q 1 ' 

whose ^-dependence is shown in Fig. 1 for q = 2 and for different values of the Levy index v. 
Smaller is the value of v faster is the convergence to D2 t k = 0, i.e. to (^-independence of 
F 2i k- This is what we expect since v — corresponds to monofractal behavior which yields 
^-independent factorial moments after a single iteration step. The v = 2 limiting case of 
the Levy index defines maximum degree of multifractality. Random cascades of this type 
are described by the log-normal law whose fixed point property discussed before is apparent 
from Fig. 1. For v > 2, which is not permitted for log-Levy laws, F^^ gets farther away 
from 5-independence. 

Our last example concerns the analogy between fractal structures and statistical mechan- 
ics [pT],pr6|| . In the thermodynamic description of multifractal cascades the rank q plays the 
role of inverse temperature fL6f . Analogously to RG calculations in statistical physics, per- 
forming the iterative sequence of transformations can be thought as varying the temperature. 
The evolution of D q ,k allows the introduction of an effective inverse temperature q^ whose 
variation with k can be obtained from the equation (q + k) v + (q — 1) k v — q(k + l) v = q^ — qk- 
The Levy index v controls the rate at which the effective inverse temperature changes in the 
course of the iteration. Therefore v plays the role of the critical exponent of the correlation 
length characterizing condensed matter systems at the critical temperature f|. 
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In conclusion, we have developed an iterative procedure capable of finding homogeneity 
rules for a distribution that describes the fluctuation properties of a certain physical system. 
A single iteration step consists of a change in the normalization point of the distribution 
compensated by a suitable change in the renormalized parameters and a rescaling. The 
number of steps required to arrive at data collapsing behavior measures how strongly the 
scaling hypothesis of type Eq. ([!]) is violated. For self-similar fluctuations it is maximal in 
the case of maximum degree of multifractality, providing the fixed point of the iteration, 
and minimal for monofractals which give rise to scaling behavior after a single iteration 
step. Concerning multiparticle dynamics, our results shed new light on how to characterize 
the degree of KNO scaling violation of the full phase-space multiplicity distributions [T?| 



in terms of novel scaling laws. Since homogeneity rules and data collapsing deserve special 
attention in many areas of physics the iterative procedure can have a number of applications 
in other fields as well. 

This research was supported by the Hungarian Science Foundation under grant OTKA- 
F4019/1992. 
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FIGURES 



v=2.1 




Iteration step k 

FIG. 1. Variation of the anomalous fractal dimension D2 with iteration step k for different 
values of the Levy index v. 
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